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Beyond the Thouless energy 
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The distribution and the correlations of the small eigenvalues of the Dirac operator are described by random 
matrix theory (RMT) up to the Thouless energy E c oc 1/s/V, where V is the physical volume. For somewhat 
larger energies, the same quantities can be described by chiral perturbation theory (chPT). For most quantities 
there is an intermediate energy regime, roughly 1/V < E < l/W, where the results of RMT and chPT agree 
with each other. We test these predictions by constructing the connected and disconnected scalar susceptibilities 
from Dirac spectra obtained in quenched SU(2) and SU(3) simulations with staggered fermions for a variety of 
lattice sizes and coupling constants. In deriving the predictions of chPT, it is important to take into account only 
those symmetries which are exactly realized on the lattice. 



The theoretical understanding of the Dirac 
eigenvalue spectrum in a finite volume has im- 
proved considerably in recent years. The small- 
est Dirac eigenvalues are described by universal 
functions which can be computed most easily in 
chiral RMT @,§]. The agreement persists up to 
the so-called Thouless energy E c which scales like 
1/L 2 , where V = L* |||-^). Beyond this energy, 
the Dirac spectrum can be described by chPT Q . 
This has been discussed in the continuum theory 
in Ref. . On a coarse lattice, the situation is dif- 
ferent, and one should take into account only the 
lattice symmetries. Here, we present an analysis 
appropriate for staggered fermions at relatively 
strong coupling and compare our predictions to 
SU(2) and SU(3) lattice gauge data. For details 
of the SU(2) analysis, we refer to Ref. ||. 

We are interested in the connected and discon- 
nected scalar susceptibilities defined by 
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respectively, where the i\k are the Dirac eigenval- 
ues and to is a valence quark mass. Most of the 
RMT-predictions for these quantities are given in 
Refs. [|],[|. The corresponding chPT-predictions 
can be derived from an effective partition func- 
tion Z by differentiating with respect to the quark 
masses ||. We consider N v generations of va- 
lence quarks of mass m v and N s generations of 
sea quarks of mass m s (corresponding to 4N V va- 
lence quarks and 4iV s sea quarks in the contin- 
uum limit). Our starting point is the following 
expression for the free energy, 

\nZ(m Vl m Sl L) oc VS(m v ,m s ) 
Q p 

where S(m v ,m s ) is the saddle-point contribu- 
tion, and the double sum represents the one- 
loop contribution coming from light composite 
bosons. The sum runs over the allowed momenta 
p (jap — lixn^jL with integer n^) and over parti- 
cle type Q (with multiplicity Kq and mass toq). 
We use p 2 = 2X^(1 _ c osp M ). 

The main task is to determine the Kq and toq 
for our particular problem. Consider first gauge 
group SU(3). The symmetry in the chiral limit is 
SV(N V + N S ) x U(l) x S\J(N V + N S ) x U(l) which 
is spontaneously broken to SU(iV„ + N s ) x U(l). 
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Figure 1. Annihilation diagram for the "flavor- 
diagonal" mesons. 

Since for staggered fermions in strong coupling 
the U(l) symmetry is anomaly-free, we expect 
(N v + N s ) 2 Goldstone bosons. The bosons made 
of different quark flavors qiQj will have a mass 
given by to 2 = A{rrn + m,j)/2. (According to the 
Gell-Mann-Oakes-Renner relation, A = 2£// 2 , 
where X = | (^)|.) We thus have N 2 —N v mesons 
of mass Am v , N 2 — N s mesons of mass Am s , and 
2N V N S mesons of mass A(m v + m s )/2. 

For "flavor-diagonal" mesons we must also con- 
sider the annihilation process in Fig. [I]. Because 
of the anomaly-free U(l) symmetry, the ampli- 
tude for qiQi — > q~jqj vanishes for m, = or 
nij = 0. Therefore, we make the following 
ansatz for the mass-squared matrix of the states 
(vivi, . . . ,V Nv V Nv ,SiSi, . . . ,s Ns s Ns ) T , 
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with an additional parameter z. The eigenvalues 
of M 2 are Am v with multiplicity N v — 1, Am s 
with multiplicity N s — 1, and X± with multiplicity 
one (the expression for X± is given in j^]). This 
completes the determination of the light boson 
spectrum of the gauge group SU(3) in Table |l|. 

For the gauge group SU(2), the symmetry in 
the chiral limit is U(2A^ + 2N S ), spontaneously 
broken to 0(2A^+2iV s ) [@. We thus have (N v + 
N S )(2N V + 2N S + 1) Goldstone particles. Some of 
the baryon {qiqj and q~iq~j) states have the same 
mass as the mesons, m 2 — A(rrii + mj)/2, giving 
rise to the light particle spectrum in Table 0. 
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Table 1 

The light particle spectrum for gauge groups 
SU(2) and SU(3). 

Table ^| determines the one-loop contribution 
to the free energy in Eq. (||). The saddle-point 
contribution is parameterized by a smooth func- 
tion of m v and m s , independent of the lat- 
tice size. Taking appropriate derivatives of In Z 
with respect to the quark masses we obtain 
the chPT-predictions for the susceptibilities of 
Eq. (|l|). In the final results, we take the limits 
m v = m s = to, N v — > 0, and N s — » 0. The fit 
parameters are A, z, and the smooth background. 
Since X can be determined independently by a fit 
to RMT, our results for the parameter A = 2E/ f 2 
also give us an estimate of 

Taking the infinite-volume limit of the chPT- 
expressions, we obtain several terms containing 
logarithms in the quark mass ||. Note, however, 
that the leading term oc In to in the chiral con- 
densate, which is expected in the quenched ap- 
proximation ||, is absent in our case because of 
the anomaly- free U(l) symmetry. 

Our results for gauge group SU(2) and SU(3) 
are displayed in Figs. || and [| The diamonds 
represent the lattice data plotted vs. the rescaled 
valence quark mass u — mVT,, the solid lines the 
(finite-volume) chPT predictions, and the dashed 
lines the RMT predictions (for topological charge 
v = 0), respectively. As expected, for u < f 2 L 2 
the data are described by RMT. For u > 1, they 
are very well described by our chPT expressions. 
(chPT breaks down for it < 1 since the p = 
modes must be treated non-perturbatively in this 
region. The deviations between lattice data and 
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Figure 2. Connected and disconnected scalar 
susceptibilities versus the rescaled valence quark 
mass for staggered fermions using gauge group 
SU(2) at (3 = A/g 2 = 2.2 and V = 8 4 and 12 4 . 



cliPT for very large u arc due to the finite lattice.) 
The domain of common applicability of RMT and 
chPT, 1 < u < f 2 L 2 , grows with the lattice size. 

In the case of the connected susceptibility in 
SU(3) (see Fig. ||) we do not see an overlap re- 
gion of RMT and chPT. The reason is that for 
this particular quantity (and also for the chiral 
condensate) the would-be leading terms both in 
RMT (for large m) and in chPT (for small m) are 
absent. This is a rather special case caused by the 
anomaly-free U(l) symmetry and by the fact that 
N s = 0. As a consequence, the Thouless energy 
for this quantity scales like 1/L 8 / 3 instead of like 
1/L 2 so that RMT breaks down for u oc L 4 / 3 . 

In conclusion, we now have a good theoreti- 
cal understanding of the finite- volume Dirac spec- 
trum also beyond the Thouless energy. Our anal- 
ysis was tailored to the case of staggered fermions 
at strong coupling where the anomaly-free U(l) 
symmetry causes the light particle spectrum to 
be different from that of the continuum theory. 

We thank M. Golterman and J.J.M. Verbaar- 
schot for helpful comments. This work was sup- 
ported in part by DFG. 
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Figure 3. Same as Fig. g but for gauge group 
SU(3) at = 6/g 2 = 5.4 and V = 6 4 and 10 4 . 
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